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In this study, 111 Cryptosporidium parvum Ild isolates from several species of animals in China, Sweden, 
and Egypt were sub typed by multilocus sequence typing (MLST). One to eleven subtypes were detected at 
each of the 12 microsatellite, minisatellite, and single nucleotide polymorphism (SNP) loci, forming 25 
MLST subtypes. Host-adaptation and significant geographical segregation were both observed in the MLST 
subtypes. A clonal population structure was seen in C parvum lid isolates from China and Sweden. Three 
ancestral lineages and the same RPGR sequence were shared by these isolates examined. Therefore, the 
present genetic observations including the higher nucleotide diversity of C. parvum lid GP60 sequences in 
Western Asia, as well as the unique distribution of lid subtypes (almost exclusively found in Asia, Europe, 
and Egypt) and in combination with the domestication history of cattle, sheep, and goats, indicated that C 
parvum lid subtypes were probably dispersed from Western Asia to other geographical regions. More 
population genetic structure studies involving various C. parvum subtype families using high-resolution 
tools are needed to better elucidate the origin and dissemination of C parvum in the world. 



Cryptosporidium parvum is a common zoonotic agent and has a wide host spectrum and geographical 
distribution, accounting for about 45% of human cryptosporidiosis cases worldwide^'^. There are extens- 
ive intraspecific variations in C. parvum populations. Currently, at least 14 C. parvum subtype families, 
including Ila, lib, lie, lld. He, llf, llg, llh. Hi, Ilk, III, Ilm, Iln, and IIo, have been identified on the basis of sequence 
analysis of the 60 KDa glykoprotein (GP60) gene^"^. Among them, Ila is the predominant subtype family in 
animals and humans worldwide whereas lid is another major zoonotic subtype family reported in Europe 
(Hungary, Germany, Portugal, Sweden, Ireland, Spain, Belgium, Romania, the United Kingdom, Netherlands, 
Slovenia, and Serbia and Montenegro), Asia (Kuwait, Iran, Jordan, India, Malaysia, and China), Egypt, and 
Australia^"^. In contrast, lie and He are anthroponotic subtype families. Other subtype families of C. parvum 
are occasionally seen in some animals or humans in the world^'^. 

The life cycle of Cryptosporidium comprises a sexual phase during which gamete differentiation and fertiliza- 
tion occurs. This process can affect the population structure of this parasite because of the possibility of genetic 
recombination between genotypes. Experimental evidence has confirmed that recombinant genotypes can 
emerge from mixed infections^. In nature, panmictic population structure representing frequent genetic exchange 
was seen in human and bovine C. parvum isolates in Scotland, Ireland, and the United States^"^\ Meanwhile, 
epidemic structure was also noticed in C. parvum human isolates in Haiti and Scotland, in which genetic exchange 
and linkage equilibrium is masked by the expansion of genetically identical parasites^^'^^. In contrast, clonal 
population structure (indicates that multilocus subtypes are relatively stable in time and place) was observed in C. 
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parvum human isolates in France^^, and bovine C. parvum in Israel 
and Italy^^'^l On the other hand, there is evidence that in one area, C. 
parvum can have an epidemic population structure in humans but a 
panmictic population structure in cattle^^, or a clonal population 
structure in humans but an epidemic population structure in ani- 
mals Therefore, the relative contribution of genetic recombination 
to the diversification of C. parvum in nature is complex, and may 
depend on diverse factors including transmission intensity, host 
movement, and host adaptation \ 

In China, C. parvum has been detected in a small number of 
animals and humans. GP60 sequence analysis showed that all C. 
parvum isolates characterized thus far belonged to Ild subtypes, 
including IldAlSGl in rodents^^, IIdA19Gl in cattle, humans, and 
urban wastewater^'^^"^^. Previously, IldAlSGl was detected in cattle, 
sheep and goats in Iran, Malaysia, and Spain^°"^^ and in humans in 
Netherlands, Australia, Iran, Malaysia, and India^'^^'^^"^^. Likewise, 
IIdA19Gl was seen in cattle in Hungary and Sweden^^'^^, and in 
humans in Portugal and Sweden^'^^. Thus, both IldAlSGl and 
IIdA19Gl are considered zonootic subtypes. 

Here we analyzed 1 1 1 C. parvum lid isolates from China, Sweden, 
and Egypt by a multilocus sequence typing (MLST) technique using 
12 microsatellite, minisatellite, and single nucleotide polymorphism 
(SNP) loci. Our findings revealed that C. parvum lid subtypes had a 
clonal population structure and probably dispersed from Western 
Asia to other geographical regions. 

Results 

MLST subtypes and sequence polymorphism. All 111 specimens 
were successfully sub typed and 1-11 sequence types were found at 
the 12 gene loci. Among them, GP60 and DZ-HRGP had relatively 
higher sequence polymorphism, with 11 and 7 subtypes being 
identified, respectively. In contrast, all 111 specimens had the same 
sequence at the normally polymorphic RPGR locus. Altogether, these 
C. parvum lid specimens formed 2S MLST subtypes (Table SI). 

Sequence data of all 12 loci were concatenated to make a multi- 
locus gene of 7209 bp in length. The genetic diversity of sequences 
was analyzed using DnaSP S. 10.01. They had 13S polymorphic 
(segregating) sites, and 12 haplotypes with a haplotype diversity 
Hd = 0.698 ± 0.039, nucleotide diversity Pi = 0.00692, and average 
number of nucleotide differences k = 47.296. The genetic diversity of 
Swedish isolates was higher than that of Chinese isolates (Table S2). 
Gene flow and genetic differentiation analysis showed significant 
differences among China, Egypt, and Sweden populations (P < 
0.001), with the pairwise fixation index (Fst) of 0.811 to 0.919. 
Thus, the value of Nm representing gene flow was only O.OS by both 
NsT and Fst estimates based on sequence data. 

Higher genetic diversity level at the GP60 gene is observed in 
Western Asia. To estimate genetic diversity level of C. parvum lid 
subtypes, the GP60 locus representing the most polymorphic 
sequence identifled thus far was used in this study. 43 GP60 sequ- 
ences belonging to IS subtypes were obtained from the GenBank 
database, including IldAHGl, IldAlSGl, IIdA16Gl, IIdA17Gl, 
IIdA18Gl, IIdA19Gl, IIdA20Gl, IIdA21Gl, IIdA22Gl, IIdA23Gl, 
IIdA24Gl, IIdA2SGl, IIdA26Gl, IIdA28Gl, and IIdA29Gl. DNA 
polymorphism was analyzed using DnaSP S. 10.01. The nucleotide 
diversity (Pi) of GP60 sequences was 0.00267 in Western Asia (Iran, 
Kuwait, and Jordan), which was higher than that in other areas, 
including European countries, other Asian countries (China, India, 
and Malaysia), Egypt, and Australia (Table S3). 

Phylogenetic analysis reveals a marked geographical separation of 
isolates. Neighbor-joining analysis of the concatenated sequences 
grouped the 2S MLST subtypes into three distinct lineages (Fig. 1). 
The lineages occupy strikingly distinct geographic distributions; 
lineage 3 occurred almost exclusively in China, with only Egyptian 
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Figure 1 | Phylogenetic relationship among MLST subtypes of C parvum 
lid isolates at 12 genetic loci as assessed by a neighbor-joining analysis of 
the concatenated nucleotide sequences, using distance calculated by the 
Kimura 2 -parameter model. The figure was generated using the softwares 
of Microsoft PowerPoint 2003 and Adobe Photoshop 8.0.1 after by the 
cycle model in MEGA version 4 opening the tree file. 

isolates (MLST-6) being included. In contrast, the Sweden isolates 
consisted of lineages 1 and 2. Likewise, the phylogenetic network 
analysis also produced the same result (Fig. SI). 

LD analysis suggests clonal population structure. To examine LD 
within and between loci, we analyzed pairwise associations of 
polymorphism across the concatenated sequences. LD was 
estimated using the ZnS statistic, which calculates the average 
correlation among alleles in all pairs of polymorphic sites. ZnS was 
elevated not only for all isolates from the three countries but also for 
isolates within China or Sweden (Table 1). These high ZnS values 
indicated the C. parvum lid subtypes had a clonal population 
structure, suggesting the absence of recombination among lineages, 
which was supported by the D' statistic (Fig. S2). In addition, the 
standardized index of association (Pa) analysis also confirmed this 
conclusion (Table S4). 

Structure analysis produces three ancestral lineages. A Bayesian 
statistical approach was used to infer population substructure in 
allelic variation using the software STRUCTURE. The 111 C. 
parvum lid isolates derived from the admixture of three ancestral 
lineages (K = 3) (Fig. S3). Cryptosporidium parvum lid isolates from 
cattle, hamsters, and horses in China and cattle in Egypt consisted of 
ancestral lineage III; whereas the 43 C. parvum lid isolates from cattle 
in Sweden formed ancestral lineages I and II (Fig. 2). Generally, 
Chinese isolates and Swedish isolates had different ancestral 
lineages. In contrast, the isolates from Egypt mainly showed the 
admixture of lineages I and III (Fig. 2). 

Discussion 

Compared to length variations of polymorphic locus by multilocus 
typing (MET) tools, multilocus sequence typing tools allow more 
definitive detection of sequence polymorphism and the inclusion 
of SNP makers\ In this study, 1-11 sequence types of C. parvum 
lid specimens were seen at each of the polymorphic loci, forming 25 
MLST subtypes (Table SI). Some ZPT, ML2, and GP60 sequences 
were restricted to a single animal species, suggesting that they maybe 
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Table 1 Linkage disequilibrium in C parvum lid isolates 

Group Sites* Polymorphic sites 


Significant 
Fisher's 


comparisons** 
X' test 


ZnS 


Rm 


Overall (China, Sweden and Egypt)^ 6836 


135 


5675/4951 


5858/5279 


0.5396 


2 


China 7031 


27 


325/0 


325/325 


0.9259 


0 


Sweden 6856 


77 


2489/1396 


2542/2349 


0.7445 


2 


*Total number of sites compared after removing gaps and missing data. 
**Number of significant comparisons without/with the Bonferroni correction. 
^LD was not calculated for Egyptian isolates due to the lack of sequence diversity. 



host adaptation^°. However, whether these genetic differences rep- 
resent different biologic characteristics is unclear. 

ZnS statistic, D' statistic, and values all confirmed that C. 
parvum lid subtypes had a clonal population structure, indicating 
that the MLST subtypes are relatively stable in time and place \ This 
population structure allows longitudinal tracking of the transmission 
of C. parvum Ild subtypes in field studies. Previously, diverse C. 
parvum structures ranging from clonal, epidemic, to panmictic have 
been reported in a few studies^^'^^. However, most of these were 
conducted by MLG analyses of polymorphism of microsatellite 
and minisatellite loci in Ila and some lie isolates. The actual genetic 
diversity of C. parvum isolates could be masked by the MLG tools, 
because many subtypes have the same length even though the 
sequences are different, such as HSP70, Mucin 1, MSC6-7, and even 
GP60 genes \ Therefore, using high-resolution MLST tools in com- 
bination with more polymorphic loci are needed before firm conclu- 
sions can be reached on the population structure of different C. 
parvum subtype families. 

Significant geographic segregation of MLST subtypes was 
observed in this study. The MLST subtypes from China, Egypt, 
and Sweden differed from each other; even though they belonged 
to the same GP60 subtype IIdA19Gl, MLST-1 and MLST-2 were 
seen in China whereas MLST-7 was noticed in Sweden (Table SI). 
Nevertheless, in Sweden, several MLST subtypes were detected in the 
same city or province, and even within the same herd (Table S5). 
Thus, even in a clonal population structure, genetically distinct C. 
parvum lid subtypes can emerge within a group of hosts. The similar 
observation was also noticed in a previous study conducted in Israel 
and Turkey^^. This may be partially attributed to the introduction of 
new animals into the herd or area via animal trade. Indeed, data 



showed that cattle trade was maintained at a considerable level in 
Sweden^ \ Therefore, frequent animal trade in Sweden would 
increase the heterogeneity of C. parvum lid subtypes on farms and 
the complexity of infections. Sequence analysis provided some evid- 
ence for the occurrence of mixed subtypes in Sweden, where a weak 
recombination event was seen by DnaSP 5.10.01 analysis (Rm = 2). 

STRUCTURE analysis revealed three ancient lineages in the C. 
parvum lid specimens examined in this study (Fig. 2), which was 
supported by the phylogenetic relationship analyses (Fig. 1 and SI). 
However, the ancestry origin of the three ancient lineages was 
unclear. In addition, all 111 specimens had the same sequence at 
the RPGR gene which had 100% of identity with another zonootic 
C. parvum subtype IlaAlSGlRl (accession no. BX538350 and 
XM_627472), whereas shared only 93.6% of homology in compar- 
ison with the anthroponotic subtype IIcA5G3 (accession no. 
EU14171 1). Nevertheless, we do not know whether RPGR is an use- 
ful marker for distinguishing zonootic and anthroponotic subtypes 
and tracking the origin of C. parvum due to the absence of genetic 
data of other subtype families. Understanding the domestication 
history of domestic animals will help us to explore the origin of C. 
parvum. Modern cattle, sheep, and goats are thought to have origi- 
nated from Western Asia (especially in the Fertile Crescent) and 
introduced to the rest of the world during the last 10,500-15,000 
years^^"^^. In addition, cattle, sheep, and goats are the most common 
hosts for C. parvum^'^. Thus, C. parvum lid family probably dispersed 
from Western Asia to other areas as a result of cattle introduction, 
gradually forming various clonal populations because of host and 
geographical segregation. However, the divergence time of different 
subtype families of C. parvum is still unclear due to the absence of 
necessary references, such as molecular clock. 
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Figure 2 | Population structure inferred by Bayesian clustering using multilocus subtype information. Ancestral population size of 3 (K = 3) was 

chosen as the best fit for the current data (Fig. 82). Upper labels are the ancestral populations; the geographic locations, MLST subtypes, and hosts of the 
isolates are at the bottom. C = cattle, H = horse, M = mouse. 
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Figure 3 | Ratio of the common C parvum subtype families found in the 
world. Numbers above bars represent numbers of samples subtyped at 
the GP60 locus. Data are based on studies pubUshed between 2004 and 
March 2013. 

In conclusion, as expected, multiple MLST subtypes of C. parvum 
lid isolates were present in various animals examined in this study. 
Host- adaptation and significant geographical segregation of MLST 
subtypes were both observed among isolates from China, Sweden, 
and Egypt. Cryptosporidium parvum Ild isolates from China, 
Sweden, and Egypt had a clonal population structure, shared the 
same RPGR sequence, and consisted of three ancestral lineages. In 
addition, the GP60 sequences of C. parvum lid isolates from Western 
Asia had a higher nucleotide diversity in comparison with that in 
other areas. Therefore, these findings, in combination with the 
unique geographical distribution of C. parvum lid subtypes (Fig. 3 
and 4) and the domestication history of cattle, sheep, and goats, 
indicated that C. parvum lid subtypes were probably dispersed from 
Western Asia to other geographical regions. Nevertheless, more 
population genetic structure studies involving various C. parvum 
subtype famiUes and high- resolution tools are needed to better elu- 
cidate the origin and spread of C. parvum in the world. 




Figure 4 | Distribution of C parvum lid subtypes in the world. The size of 
black circle indicates the number of isolates. Blue and red circles represent 
isolates from animals and humans, respectively. The figure was generated 
using the softwares of Maplnfo Professional 7.0, Microsoft PowerPoint 
2003 and Adobe Photoshop 8.0.1. 



Methods 

Cryptosporidium specimens. A total of 1 1 1 C. parvum lid specimens were used in 
this study including Chinese isolates from dairy cattle, horses, and hamsters and 
Swedish and Egyptian isolates from dairy cattle (Table SI). Most of the specimens 
were identified in previous studies*'^*^'^*^, whereas the 10 bovine isolates from Egypt and 
six horse and five bovine isolates from China were diagnosed as C. parvum lid 
subtypes by DNA sequence analysis of the GP60 gene^^. 

Multilocus sequence typing. Twelve micro satellite, minisatellite, and SNP loci were 
used in the present study. Primers and amplification conditions used in nested- PCR 
analysis of MSC6-7 (serine repeat antigen), RPGR (retinitis pigmentosa GTPase 
regulator), MSC6-5 (hypothetical trans -membrane protein), DZ-HRGP 
(hydroxyproline-rich glycoprotein), ZPT (zinc protease telomerase), CP56 (56 kDa 
trans -membrane protein), Chom3T (T-rich gene fragment), CP47 (47 kDa protein), 
HSP70 (70 kDa heat shock protein). Mucin 1 (mucin-like protein), ML2 (possible 
sporozoite cystein-protein), and GP60 genes were previously described^ KOD-FX- 
Neo amplification enzyme (Toyobo Co. Ltd, Osaka, Japan) was used in PCR. 400 ng/ 
|il of non-acetylated bovine serum albumin (Solarbio Co. Ltd, Beijing, China) was 
used in primary PCR to neutralize PCR inhibitors. The secondary PCR products were 
examined by agarose gel electrophoresis and visualized after GelRed™ (Biotium Inc., 
Hayward, CA) staining. The secondary PCR products were sequenced on an ABI 
3730 DNA Analyzer (Applied Biosystems, Foster City, USA), using the secondary 
primers and the Big Dye Terminator v3. 1 Cycle Sequencing kit (Applied Biosystems). 
The sequence accuracy was confirmed by bi-directional sequencing and by 
sequencing a new PCR product if necessary. 

Phylogenetic relationship and DNA polymorphism analysis. Sequence alignment 
was made using the program ClustalX 1.83 (http://www.clustal.org/). Neighbor 
joining trees were constructed using the program Phylip 3.69, based on the 
evolutionary distances calculated by Kimura-2-parameter model. The reliability of 
these trees was assessed using the bootstrap analysis with 1,000 replicates. The 
phylogenetic network analysis was done using the median -joining algorithm as 
implemented in Network 4.611 based on the polymorphism of multilocus sequence 
data from 1 1 1 C. parvum lid isolates^^ DNA polymorphism, pairwise Est values for 
the different populations, and gene flow were calculated using DnaSP 5. 10. OP''. 

Linkage disequilibrium. Sequences from 12 loci were combined in a single contig 
and analyzed for linkage disequilibrium (LD) across the entire composite sequence. 
LD analysis was conducted using software DnaSP 5.10.0P''. D' was calculated for all 
pairs of sites. Both the two-tailed Eisher's exact test and the test were used to 
determine significance of the associations between polymorphic sites. The average LD 
was estimated by using the ZnS statistic, which averages LD over all pairwise 
comparisons for S polymorphisms in N sequences. Alternatively, as a confirmation 
for the ZnS statistic, the standardized index of association (/^a) by Habould and 
Hudson (2000) was also used to assess the LD across all loci. The index and its 
probability under a null model of complete panmixia were calculated using LIAN 
version 3.5'*° with hypothesis testing by a parametric method. 

Structure analysis. The subpopulations of C. parvum lid subtypes were determined 
based on the allelic data using STRUCTURE 2.3.3 by K-means partitional clustering 
and the admixture modeP\ Twenty simulation runs were conducted for each of K = 2 
to 10 using a length of burn-in of 10^ and 10^ replicates of Markov chain Monte Carlo 
simulation. The average membership coefficients for the 20 simulation runs of a given 
K value were generated by CLUMPP vl.1.2 (windows version)*^ and a graphical 
presentation of the average membership coefficients for each isolate was generated in 
Microsoft Excel. The most appropriate number of K was calculated using an ad hoc 
statistic-based approach implemented in software Structure Harvester vO.6.93, as 
described previously'^^ 

Nucleotide sequence accession numbers. Representative nucleotide sequences were 
deposited in the GenBank under accession numbers KC823118 to KC823157. 



1. Xiao, L. & Ryan, U. [Molecular epidemiology] Cryptosporidium And 
Cryptosporidiosis [Layer, R. &Xiao, L. (eds)] [119-150] (CRC Press and IWA 
Publishing, Boca Raton, 2008). 

2. Xiao, L. Molecular epidemiology of cryptosporidiosis: an update. Exp. Parasitol. 
124, 80-89 (2010). 

3. Insulander, M. et al. Molecular epidemiology and clinical manifestations of 
human cryptosporidiosis in Sweden. Epidemiol. Infect. 141, 1009-1020 (2013). 

4. Wang, R. J. et al. Characteristics of Cryptosporidium transmission in preweaned 
dairy cattle in henan, China. /. Clin. Microbiol. 49, 1077-1082 (2011). 

5. Amer, S. et al. Identity and public health potential of Cryptosporidium spp. in 
water buffalo calves in Egypt. Vet. Parasitol. 191, 123-127 (2012). 

6. Iqbal, A., Lim, Y. A., Surin, J. & Sim, B. L. High diversity of Cryptosporidium 
subgenotypes identified in Malaysian HIV/ AIDS individuals targeting gp60 gene. 
PLoS One 7, e31139 (2012). 

7. Imre, K. et al. Zoonotic Cryptosporidium parvum in Romanian newborn lambs 
(Ov/s aries). Vet. Parasitol. 191, 119-122 (2013). 

8. Eeng, X., Rich, S. M., Tzipori, S. & Widmer, G. Experimental evidence for genetic 
recombination in the opportunistic pathogen Cryptosporidium parvum. Mol. 
Biochem. Parasitol. 119, 55-62 (2002). 



SCIENTIFIC REPORT! | 4 : 4208 | DOI: 1 0.1 038/srep04208 



4 



9. Mallon, M. et al. Population structures and the role of genetic exchange in the 
zoonotic pathogen Cryptosporidium parvum. J. Mol. Evol. 56, 407-417 (2003). 

10. Herges, G. R. et al. Evidence that Cryptosporidium parvum populations are 
panmictic and unstructured in the Upper Midwest of the United States. Appl. 
Environ. Microbiol. 78, 8096-8101 (2012). 

1 1 . De Waele, V. et al. Panmictic structure of the Cryptosporidium parvum population 
in Irish calves: influence of prevalence and host movement. Appl. Environ. 
Microbiol. 79, 2534-2541 (2013). 

12. Mallon, M. E., Macleod, A., Wastling, J. M., Smith, H. & Tait, A. Multilocus 
genotyping of Cryptosporidium parvum Type 2, population genetics and sub- 
structuring. Infect. Genet. Evol. 3, 207-218 (2003). 

13. Ngouanesavanh, T. et al. Cryptosporidium population genetics: evidence of 
clonality in isolates from France and Haiti. /. Eukaryot. Microbiol. 53(Suppl), 
33-36 (2006). 

14. Tanriverdi, S. & Widmer, G. Differential evolution of repetitive sequences in 
Cryptosporidium parvum and Cryptosporidium hominis. Infect. Genet. Evol. 6, 
113-122 (2006). 

15. Drumo, R. et al. Evidence of host-associated populations of Cryptosporidium 
parvum in Italy. Appl. Environ. Microbiol. 78, 3523-3529 (2012). 

16. Lv, C. C. et al. Cryptosporidium spp. in wild, laboratory, and pet rodents in china: 
prevalence and molecular characterization. Appl. Environ. Microbiol. 75, 
7692-7699 (2009). 

17. Li, N. et al. Molecular surveillance of Cryptosporidium spp., Giardia duodenalis, 
and Enterocytozoon bieneusi by genotyping and subtyping parasites in 
wastewater. PLoS Negl. Trop. Dis. 6, el809 (2012). 

18. Zhang, W. Z. et al. Distribution and genetic characterizations of Cryptosporidium 
spp. in pre-weaned dairy calves in Northeastern China's Heilongjiang Province. 
PLoS One 8, e54857 (2013). 

19. Wang, L. et al. Zoonotic Cryptosporidium species and Enterocytozoon bieneusi 
genotypes in HIV-positive patients on antiretroviral therapy. /. Clin. Microbiol. 
51, 557-563 (2013). 

20. Quilez, J. et al. Cryptosporidium species and subtype analysis from dairy calves in 
Spain. Parasitology 135, 1613-1620 (2008). 

21. Quilez, J. et al. Cryptosporidium genotypes and subtypes in lambs and goat kids in 
Spain. Appl. Environ. Microbiol. 74, 6026-6031 (2008). 

22. Nazemalhosseini-Mojarad, E. et al. Subtype analysis of Cryptosporidium parvum 
and Cryptosporidium hominis isolates from humans and cattle in Iran. Vet. 
Parasitol. 179, 250-252 (2011). 

23. Muhid, A., Robertson, I., Ng, J. & Ryan, U. Prevalence of and management factors 
contributing to Cryptosporidium sp. infection in pre-weaned and post-weaned 
calves in Johor, Malaysia. Exp. Parasitol. 127, 534-538 (2011). 

24. Wielinga, P. R. et al. Molecular epidemiology of Cryptosporidium in humans and 
cattle in The Netherlands. Int. J. Parasitol. 38, 809-817 (2008). 

25. Ng, J., MacKenzie, B. & Ryan, U. Longitudinal multi-locus molecular 
characterisation of sporadic Australian human clinical cases of cryptosporidiosis 
from 2005 to 2008. Exp. Parasitol. 125, 348-356 (2010). 

26. Ajjampur, S. S. et al. Multisite study of cryptosporidiosis in children with diarrhea 
in India. /. Clin. Microbiol. 48, 2075-2081 (2010). 

27. Plutzer, J. & Karanis, P. Genotype and subtype analyses of Cryptosporidium 
isolates from cattle in Hungary. Vet. Parasitol. 146, 357-362 (2007). 

28. Silverlas, C., Bosaeus-Reineck, H., Naslund, K. & Bjorkman, C. Is there a need for 
improved Cryptosporidium diagnostics in Swedish calves? Int. J. Parasitol. 43, 
155-161 (2013). 

29. Alves, M., Xiao, L., Antunes, F. & Matos, O. Distribution of Cryptosporidium 
subtypes in humans and domestic and wild ruminants in Portugal. Parasitol. Res. 
99, 287-292 (2006). 

30. Xiao, L. etal. Host adaptation and host-parasite co-evolution in Cryptosporidium: 
implications for taxonomy and public health. Int. J. Parasitol. 32, 1773-1785 
(2002). 

31. Widgren, S. & Frossling, J. Spatio-temporal evaluation of cattle trade in Sweden: 
description of a grid network visualization technique. Geospat. Health 5, 1 19-130 
(2010). 



32. Gotherstrom, A. et al. Cattle domestication in the Near East was followed by 
hybridization with aurochs bulls in Europe. Proc. Biol. Sci. 272, 2345-2350 (2005). 

33. Beja-Pereira, A. et al. The origin of European cattle: evidence from modern and 
ancient DNA. Proc. Natl. Acad. Sci. U. S. A. 103, 8113-8118 (2006). 

34. Fernandez, H. et al. Divergent mtDNA lineages of goats in an Early Neolithic site, 
far from the initial domestication areas. Proc. Natl. Acad. Sci. U. S. A. 103, 
15375-15379 (2006). 

35. Chessa, B. et al. Revealing the history of sheep domestication using retrovirus 
integrations. Science 324, 532-536 (2009). 

36. Ajmone-Marsan, P., Garcia, J. F. & Lenstra, J. A. On the origin of cattle: How 
aurochs became cattle and colonized the world. Evolutionary Anthropology: 
Issues, News, and Reviews 19, 148-157 (2010). 

37. Sulaiman, I. M. et al. Unique endemicity of cryptosporidiosis in children in 
Kuwait. /. Clin. Microbiol. 43, 2805-2809 (2005). 

38. Bandelt, H. J., Forster, P. & Rohl, A. Median-joining networks for inferring 
intraspecific phylogenies. Mol. Biol. Evol. 16, 37-48 (1999). 

39. Rozas, J., Sanchez-DelBarrio, J. C, Messeguer, X. & Rozas, R. DnaSP, DNA 
polymorphism analyses by the coalescent and other methods. Bioinformatics 19, 
2496-2497 (2003). 

40. Haubol, B. & Hudson, R. R. LIAN 3.0: detecting linkage disequilibrium in 
multilocus data. Bioinformatics 16, 847-848 (2000). 

41. Falush, D., Stephens, M. & Pritchard, J. K. Inference of population structure using 
multilocus genotype data: Linked loci and correlated allele frequencies. Genetics 
164, 1567-1587 (2003). 

42. Jakobsson, M. & Rosenberg, N. A. CLUMPP: A cluster matching and permutation 
program for dealing with label switching and multimodality in analysis of 
population structure. Bioinformatics 23, 1801-1806 (2007). 

43. Earl, D. A. & vonHoldt, B. M. STRUCTURE HARVESTER: a website and 
program for visualizing STRUCTURE output and implementing the Evanno 
method. Conservation Genet. Resour. 4, 359-361 (2012). 



Acknowledgments 

This study was supported in part by the State Key Program of National Natural Science 
Foundation of China (31330079), the National Natural Science Foundation of China 
(U1204328, 31172311, 31302079), the International Cooperation and Exchange Funds of 
the National Natural Science Foundation of China (31110103901), and the Key National 
Science and Technology Specific Projects (2012ZX10004-220). 

Author contributions 

Conceived and designed the experiments: L.X.Z., L.H.X. and R.J.W. Performed the 
experiments: R.J.W., F.C.J., M.Q. and H.J.D. Analyzed the data: R.J.W., G.Q.L., C.C.L., 
Z.F.Z., H.L.W. and Y.R.S. Contributed reagents/materials/analysis tools: C.A., C.B., S.A., 
A.Q.L., Y.Y.F. and C.S.N. Wrote the paper: R.J.W., L.X.Z. and L.H.X. 

Additional information 

Supplementary information accompanies this paper at http://www.nature.com/ 
scientificreports 

Competing financial interests: The authors declare no competing financial interests. 

How to cite this article: Wang, R.J. et al. Cryptosporidium parvum lid family: clonal 
population and dispersal from Western Asia to other geographical regions. Sci. Rep. 4, 4208; 
DOL10.1038/srep04208 (2014). 



©0®© 



This work is licensed under a Creative Commons Attribution- 
NonCommercial-NoDerivs 3.0 Unported license. To view a copy of this license, 



visit http://creativecommons.Org/licenses/by-nc-nd/3.0 



SCIENTIFIC REPORTS | 4 : 4208 | DOI: 1 0.1 038/srep04208 



5 



